Monopole Quantum Numbers in the Staggered Flux Spin Liquid 
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Algebraic spin liquids, which are exotic gapless spin states preserving all microscopic symmetries, have been 
widely studied due to potential realizations in frustrated quantum magnets and the cuprates. At low energies, 
such putative phases are described by quantum electrodynamics in 2 + 1 dimensions. While significant progress 
has been made in understanding this nontrivial interacting field theory and the associated spin physics, one 
important issue which has proved elusive is the quantum numbers carried by so-called monopole operators. 
Here we address this issue in the "staggered-flux" spin liquid which may be relevant to the pseudogap regime 
in high-Tc. Employing general analytical arguments supported by simple numerics, we argue that proximate 
phases encoded in the monopole operators include the familiar Neel and valence bond solid orders, as well 
as other symmetry-breaking orders closely related to those previously explored in the monopole-free sector of 
the theory. Surprisingly, we also find that one monopole operator carries trivial quantum numbers, and briefly 
discuss its possible implications. 
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I. INTRODUCTION 

When frustration or doping drives quantum fluctuations 
sufficiently strong to destroy symmetry-breaking order even 
at zero temperature, exotic ground states known as spin liq- 
uids emerge. "Algebraic spin liquids" comprise one class 
in which the spins appear "critical", exhibiting gapless ex- 
citations and power-law correlations which, remarkably, can 
be unified for symmetry-unrelated observables such as mag- 
netic and valence bond solid fluctuations. This unification of 
naively unrelated correlations is a particularly intriguing fea- 
ture, in part because it constitutes a "smoking gun" prediction 
for the detection of such phases. 

While the unambiguous experimental observation of a 
quantum spin liquid (either gapless, or the related topological 
variety) remains to be fulfilled, there are a number of candi- 
date materials which may host such exotic ground states. Re- 
cently the spin- 1/2 kagome antiferromagnet herbertsmithite 
has emerged as a prominent example ^ 1 ' 2 ' 3 ' 4 ' 5,6 and several gap- 
less spin liquid proposals 7,8,9 ' 10,11 , as well as a more conven- 
tional valence bond solid phas e - 12,13,14 ' 15 , have been put forth 
for this material. Furthermore, the cuprates have long been 
speculated to harbor physics connected to an algebraic spin 
liquid — the so-called "staggered-flux" state which we will fo- 
cus on here — in the pseudogap regime of the phase diagram 
(for a recent comprehensive review, see Ref.llq). 

On the theoretical end, our understanding of algebraic 
spin liquids has grown dramatically over the past several 
years. Such states are conventionally formulated in terms 
of fermionic, charge-neutral "spinon" fields coupled to a 
U(l) gauge field, whose low-energy dynamics is described 
by compact quantum electrodynamics in 2 + 1 dimensions 
(QED3). Much effort has been focused on addressing two 
basic questions concerning these states. First, can they be 
stable? In more formal terms, is criticality in QED3 pro- 
tected, or are there relevant perturbations allowed by symme- 
try which genetically drive the system away from the critical 
fixed point? And second, if algebraic spin liquids are stable, 
what are the measurable consequences for the spin system? 



Both are nontrivial questions that require consideration of two 
classes of operators in QED3 — those that conserve gauge flux 
such as spinon bilinears, and "monopole operators" that incre- 
ment the gauge flux by discrete units of 2-7T. 

While QED3 is known to be a strongly interacting field the- 
ory which lacks a free quasi-particle description, the theory 
can nevertheless be controlled by generalizing to a large num- 
ber N of spinon fields and performing an analysis in powers of 
1/N. Within such a large- A" approach, the answer to the first 
question has been rigorously shown to be 'yes' — such phases 
can in principle be stable. I7 In particular, despite some contro- 
versy concerning the relevance of monopoles, it has now been 
established that such operators are strongly irrelevant in the 
large- N limit, their scaling dimension scaling linearly with 

Significant progress has also been made in addressing the 
second question, particularly in the monopole-free sector. The 
effective low-energy QED3 theory for algebraic spin liquids 
is known to possess much higher symmetry than that of the 
underlying microscopic spin Hamiltonian, leading to the re- 
markable unification of naively unrelated competing orders 
noted above. Furthermore, the machinery of the projective 
symmetry group 19 allows one to establish how correlations of 
flux-conserving operators in QED3 relate to physical observ- 
ables such as Neel or valence bond solid correlations — , and 
the large- N analysis additionally provides quantitative predic- 
tions for the corresponding scaling dimensions 21 . 

The physical content of monopole operators in QED3, how- 
ever, is much less understood. Essentially, the difficulty here 
is that, due to gauge-invariance, determining monopole quan- 
tum numbers requires examination of full many-body spinon 
wavefunctions, rather than just a few low-energy single- 
particle states as suffices, say, for the spinon bilinears. Al- 
though the monopoles are highly irrelevant in the large- N 
limit, their scaling dimensions may become of order unity for 
realistic values of N (e.g., N = 4 for the staggered-flux state), 
so understanding the competing orders encoded in these op- 
erators becomes an important and physically relevant issue. 
Moreover, since monopoles are allowed perturbations in com- 



pact QED3 which can in principle destroy criticality for small 
enough N, one would like to identify the leading symmetry- 
allowed monopole operators. Some progress on these issues 
has been made for gapless spin liquids on the triangular and 
kagome lattices 8,10 ' 22 ' 23 ' 24 , though in the important staggered- 
flux state the physics encoded in the monopoles remains a 
mystery 20,25 . The goal of this paper is to generalize the tech- 
niques employed earlier in the former cases to deduce the 
monopole quantum numbers for the staggered-fiux state and 
reveal the competing orders encoded in this sector of the the- 
ory. 



theory on the square lattice, 
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where c r -m are spinful fermionic operators, the first sum in 
Eq. (0 represents a lattice curl summed over all plaquettes, 
and the divergence of the electric field E rr i is constrained such 
that 



A. Assumptions and Strategy 
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C rrv C r 
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Let us at the outset discuss the core assumptions on which 
our quantum-number analysis will be based. First, we will 
assume that it is sufficient to study monopoles at the mean- 
field level. That is, we will treat the flux added by a monopole 
operator as a static background "felt" by the spinons. This is 
reasonable coming from the large-iV limit, where gauge fluc- 
tuations are strongly suppressed, and is in fact the standard ap- 
proach adopted when discussing such flux insertions (see, e.g., 
Ref. [la) . The second, and more crucial, assumption we em- 
ploy is that the quantum numbers for the leading monopoles 
(those with the slowest-decaying correlations) can be obtained 
from the difference in quantum numbers between the mean- 
field ground states with and without the flux insertion. Put 
more physically, the leading monopole quantum numbers are 
taken to be the momentum, angular momentum, etc. , imparted 
to the spinon ground states upon flux insertion. 

The latter is equivalent to assuming that 1.) the flux in- 
sertion is "adiabatic" in the sense that the fermionic spinons 
remain in their relative ground state everywhere between the 
initial and final state and 2.) no Berry phases are accumu- 
lated during this evolution. The first point follows because 
if the fermions remain in their relative ground before and af- 
ter the flux insertion, then this ought to be true everywhere in 
between as well. Such an assumption is quite delicate given 
that the mean-field states we will study are gapless in the ther- 
modynamic limit. We will not attempt to justify this point 
rigorously, but we note that treating the problem in this way is 
in the same spirit as the conventional mean-field treatment of 
flux insertions mentioned above. If invalid, then treating flux 
insertions as a static background in the first place may not be 
a very useful starting point for addressing this problem. 

Assuming no Berry phases is equally delicate. It is worth 
mentioning that this assumption is known to break down in 
certain cases. As an illustration, consider the following gauge 



The standard electric-magnetic duality can be applied in the 
limit v/t — > oor& in which case one obtains a pure gauge the- 
ory with (A • E) r = (—l) rm+ry . Such an analysis reveals 
that the leading monopole operators carry nontrivial quan- 
tum numbers as a consequence of Berry phase effects, 26 even 
though the quantum numbers of the fermions clearly can not 
change in this limit. The root of these nontrivial quantum 
numbers can be traced to the fact that the electric field diver- 
gence changes sign between neighboring sites. If one alterna- 
tively considered a pure gauge theory with vanishing electric 
field divergence, then no such Berry phases arise. Since in 
the staggered-fiux state of interest the physical Hilbert space 
has exactly one fermion per site and thus a vanishing electric 
field divergence, we believe that it is reasonable to suspect that 
Berry phases do not play a role there as well. 

Given these assumptions, we will adopt the following strat- 
egy below. First, we will give a quick overview of the x- 
flux and staggered flux states, deriving a low-energy mean- 
field Hamiltonian for these states as well as the symmetry 
properties for the continuum fields. We will then consider 
±27r flux insertions, and in particular obtain the transforma- 
tion properties for the four quasi-localized zero-modes which 
appear. Armed with this information, we will follow closely 
the monopole study of Refs. |2j and |23| and constrain the 
monopole quantum numbers as much as possible using vari- 
ous symmetry relations which must generically hold on phys- 
ical states, such as two reflections yielding the identity. The 
ambiguities that remain will be sorted out by appealing to 
general quantum number conservation and simple numerical 
diagonalization for systems with convenient geometries and 
gauge choices. This will allow us to unambiguously deter- 
mine the monopole quantum numbers, subject to the above 
assumptions. We will then explore the competing orders en- 
coded in the monopole operators, and close with a brief dis- 
cussion of some outstanding questions. 



II. PRELIMINARIES 

A. Overview of 7r-flux and staggered-flux states 

Although we will ultimately be interested in exploring 
monopole quantum numbers in the staggered-flux state, we 
will use proximity to the 7r-flux state in our analysis and thus 
discuss both states here. Consider, then, a square-lattice anti- 
ferromagnet with Hamiltonian 
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Mean field descriptions of the 7r-flux and staggered-flux states 
can be obtained from Eq. (O by first decomposing the spin 
operators in terms of slave fermions via 



S r - ^fla^apfrp, 



(6) 



where cr is a vector of Pauli spin matrices and the fermions 
are constrained such that there is exactly one per site. As dis- 
cussed in Refs. [l9ll27ll28l there is an SU(2) gauge redundancy 
in this rewriting. The resulting bi-quadratic fermion Hamil- 
tonian can then be decoupled using a Hubbard-Stratonovich 
transformation, giving rise to a simple free-fermion Hamilto- 
nian at the mean-field level of the form 
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The 7r-flux state corresponds to an ansatz where the 
fermions hop in a background of tt flux per plaquette; i.e., a rr / 
is chosen such that (Ax a) = tt around each square. This state 
retains the full SU(2) gauge redundancy inherent in Eq. (0. 
As the name suggests, the staggered-flux state corresponds to 
an ansatz in which the fermions hop in flux which alternates in 
sign between adjacent plaquettes; i.e., (A x a) = ±$, where 
$ is the flux magnitude. Note that this ansatz reduces to the 
7r-flux ansatz when <J> = tt since tt flux and —tt flux are equiv- 
alent on the lattice. In contrast to the 7r-flux state, there is only 
a U(l) gauge redundancy remaining here. Note also that, de- 
spite appearances, staggering the flux does not break transla- 
tion symmetry. Rather, translation symmetry (and others) are 
realized nontrivially as a result of gauge redundancy — the op- 
erators transform under a projective symmetry group 19 . Both 
ansatzes in fact preserve all microscopic symmetries of the 
original spin Hamiltonian, namely, x and y translations T x>y , 
ir/2 rotations about plaquette centers R v / 2 , x-re flection about 
square lattice sites R x , time reversal T, and SU(2) spin sym- 
metry. Notably, there is no symmetry leading to conservation 
of gauge flux, which is why monopole operators are in princi- 
ple allowed perturbations. 



B. Continuum Hamiltonian and symmetry transformations 

To derive a continuum Hamiltonian and deduce how the 
fields transform under the microscopic symmetries, we will 



now choose a gauge and set e Mrr ' = 1 on vertical links and 
e m rr' — (—l) y on horizontal links. Although this corresponds 
to tt flux, the transformations for the staggered-flux state can 
still be readily obtained from this choice. Furthermore, adopt- 
ing this starting point yields the same continuum Hamiltonian 
as if we had chosen a staggered-flux pattern, up to irrelevant 
perturbations. 20 To obtain the spectrum we take a two-site unit 
cell and label unit cells by vectors R = n x St + 2n y y (n x . y are 
integers) which point to sites on sublattice 1; sublattice 2 is 
located at R + y. We denote the spinon operators on the two 
sublattices by /r q i,2, where a labels spin. The band structure 
is straightforward to evaluate, and at the Fermi level one finds 
two Dirac points at momenta ±Q, with Q = (tt/2, tt/2). Fo- 
cusing on low-energy excitations in the vicinity of these Dirac 
points, a continuum theory can be derived by expanding the 
lattice fermion operators as follows, 
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Here we have introduced four flavors of two-component Dirac 
fermions 4>aA, where a labels the spin and A = R/L labels 
the node. We then obtain the continuum mean-field Hamilto- 
nian 
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MF 
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(10) 



where v ~ t is the Fermi velocity and t% are Pauli matrices 
that contract with the Dirac indices. 

It is a straightforward exercise to deduce the transformation 
properties of continuum fields from Eqs. ([8]) and {9). For ei- 
ther the 7r-fiux or staggered-flux states, these can be realized 
as follows: 
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where in addition to the spin and Dirac matrices we have intro- 
duced Pauli matrices n\ B that contract with the node indices. 
In the 7r-fiux and staggered-flux cases, these transformations 
can be followed by an arbitrary SU(2) and U(l) gauge trans- 
formation, respectively. For the former, it will prove useful to 
consider a particle-hole gauge transformation Cq which is an 
element of the SU(2) gauge group and transforms the lattice 
fermion operators as 
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It follows that for the continuum fields we have 
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We stress that in the staggered-flux state Cg reverses the sign 
of the flux microscopically and therefore does not represent a 
valid gauge transformation there. 



C. Flux insertion and zero-modes 

Next we discuss the sector of the theory with ±2ir flux in- 
serted over a large area compared to the lattice unit cell. Treat- 
ing the flux as a static background, the mean-field Hamilto- 
nian then becomes 



H 



MF,q 



iv^[(d x - ial)r x + {d y - ic&)T»W. (19) 



The vector potential is chosen such that V x a q = 2irq, where 
q = ±1 is the monopole charge. It is well known that the 
above Hamiltonian admits one quasi-localized zero-mode for 
each fermion flavor, 29 four in this case. These zero-modes can 
be obtained by replacing ^ Qj 4(x) — v 4> a A,q('x)d a A,q, where 
<\>aA,q (x) is the quasi-localized wavefunction and d a A,q an- 
nihilates the corresponding state. Employing the Coulomb 
gauge, the wave functions are simply 
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It follows that the zero-mode operators d a A,q transform in ex- 
actly the same way as ip a Aj, so the transformations can be 
read off from Eqs. (fTTT i through (fT~5T > and ( fT~8b . For example, 
under reflections, we have d a R/L,q ~* iqd a L/R,-q- 

Since gauge-invariant states are half-filled, two of the four 
zero-modes must be filled in the ground states here. Thus, it 
will be convenient to introduce the following short-hand nota- 
tion: 
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d\R,qdiR, q + dlL,qdlL,q (22) 
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Of these, £>i,2,3 are spin-singlets, while 1)4,5,6 are s P m 
triplets. The transformation properties of these operators un- 
der the microscopic symmetries, as well as the gauge trans- 
formation Cg in the case of the 7r-flux state, are given in Table 
U Note that Cg changes the sign of the monopole charge q, 
indicating that the states with +2ir flux and — 2tt flux are not 
physically distinct in the 7r-flux case. We will use this fact to 
infer which of the leading monopole operators have dominant 
amplitudes in the neighboring staggered-flux state in Sec.lIVI 
We pause now to comment in greater detail on the subtlety 
with determining the staggered-flux monopole quantum num- 
bers. Naively, one might suspect that these can be inferred 
from the transformation properties of the zero-modes, which 



TABLE I: Transformation properties of the operators Dj tq defined 
in Eqs. J22b through d27l l which fill two of the four zero-modes in 
the presence of a 2irq flux insertion. The gauge transformation Co 
applies only in the 7r-flux state. 
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we have at hand. Realizing the microscopic symmetries, how- 
ever, generically requires gauge transformations, which leads 
to inherent ambiguities in how the fields transform. In partic- 
ular, for the staggered-flux case, there is an arbitrary overall 
U(l) phase in the transformations quoted in Table U and a still 
greater ambiguity in the 7r-flux state due to its larger SU(2) 
gauge group. But the monopole operators are gauge-invariant, 
so one must instead examine the symmetries of the full many- 
body wavefunctions, which are gauge invariant, rather than 
single-particle states. In what follows we will first deduce 
the transformation properties of flux insertion operators $t 
which add 2ixq flux to the ground state and fill two of the zero 
modes, 
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Here \q) represents the filled Dirac sea in the presence of 2ixq 
flux with all four zero-modes empty and |0) is the ground state 
in the absence of a flux insertion. The monopoles we will ulti- 
mately be interested in will be simply related to these objects. 
Once we know the transformation properties of $^ it will be 
trivial to read off the monopole quantum numbers. 



III. QUANTUM NUMBER DETERMINATION 
A. Symmetry relations 

As a first step, we will now constrain the quantum num- 
bers of the operators $J defined above using various symme- 
try relations which must hold when acting on gauge-invariant 
states. In particular, we will utilize the following, 
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Furthermore, all lattice symmetries must commute with time- 
reversal (when acting on gauge-invariant states). 

Quite generally, we expect the following transformations to 



hold, 
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The last transformation holds only for the 7r-fiux state. All 
phases introduced above are arbitrary at this point, but will 
be constrained once we impose symmetry relations on gauge- 
invariant states which have two of the zero-modes filled. 
Moreover, since time-reversal is anti-unitary, we have chosen 
the phases of \q) such that no additional phase factor appears 
under this symmetry. 

Consider reflections first. Equation d29l ) and commuta- 
tion with time-reversal imply that e 10 * = s, for some q- 
independent sign s. The value of s is insignificant, however, 
since we can always remove it by sending |+) — ► s\+). Hence 
we will take 
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For translations, Eqs. (f30b and PIT ), as well as commutation 
with time-reversal, yield 
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for some unknown signs s x ,y Similarly, Eq. (1321 and com- 
mutation with time-reversal allow us to determine 9 9 /2 up to 

signs sl /2 : 
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Let us turn now to the 7r-flux state, where the mean-field 
Hamiltonian is invariant under the particle-hole transforma- 
tion Cg as well. For the moment we will treat this operation 
like the other physical symmetries, which is merely a conve- 
nient trick for backing out the quantum numbers of interest 
for the staggered-fiux state. In particular, we will assert that 
C G = 1 and that this particle-hole transformation commutes 
with the physical symmetries when acting on half-filled states. 
This yields 



SG, 



(42) 



for an undetermined sign sg, and also gives the useful con- 
straint 



S n/2 — S ?r/2- 



(43) 



It follows from the last equation that 



Since the staggered-flux mean-field continuously connects to 
the 7r-flux ansatz, we will assume that the latter two con- 
straints hold in the staggered-flux case as well. (We could al- 
ternatively obtain this result using the numerics from the next 
section, without appealing to the 7r-flux state.) 

To recap, in our study of the flux-insertion operators $j 
thus far, we have shown that symmetry relations highly con- 
strain how these objects transform, and proximity to the x- 
flux state constrained these transformations even further. All 
that remains to be determined are the signs s x and s + /2 which 
appear under ^-translations and tt/2 rotations. In the follow- 
ing section we argue that these can be obtained by employing 
general quantum number conservation arguments supported 
by simple numerical diagonalization. 



B. Numerical Diagonalization 

To determine the remaining signs s x and s„/ 2 , we will now 
discuss our numerical diagonalization study of the mean-field 
Hamiltonian with and without a flux insertion, and discuss a 
more intuitive quantum number conservation argument which 
is consistent with these numerics. The basic idea behind our 
numerics is that we will judiciously choose the system geom- 
etry and gauge such that the symmetry under consideration 
can be realized without implementing a gauge transformation. 
This is a crucial point, as only in this case can we avoid over- 
all phase ambiguities that would otherwise appear in such a 
mean-field treatment. Once the single-particle wavefunctions 
with and without a flux insertion are at hand, one can proceed 
to deduce the transformation properties of the corresponding 
many-body wavefunctions and, in turn, the flux-insertion op- 
erators <1>J by using the results of the previous section. 

Consider first tt/2 rotations. Here we diagonalize the mean- 
field Hamiltonian in a square L by L system with open bound- 
ary conditions and L odd so that the system is invariant under 
tt/2 rotations about the central plaquette's midpoint. For all 
flux configurations we choose a rotationally symmetric gauge 
so that tt/2 rotations are realized trivially. We work in the 
7r-flux ansatz for simplicity, though staggering the flux can 
easily be done and clearly does not change any of the results. 
Flux is inserted over the few innermost "rings" of the sys- 
tem, and the "zero-modes" that appear quasi-localized around 
the 2tt flux can be unambiguously identified by examining the 
spread of their wave functions. (The "zero-modes" here are 
pushed away from zero energy due to finite-size effects; for 
each spin, one is pushed to higher energy while the other to 
lower energy.) We consider a variety of system sizes, with up 
to on the order of 1000 lattice sites, and obtain consistent re- 
sults in all cases examined. (More details on these numerics 
can be found in Ref . |2J, which carried out a similar study on 
the triangular lattice.) 

In particular, by considering the six ways of filling the zero- 
modes, we find numerically that there are four —1 and two 
+1 rotation eigenvalues for the operators $t To then back 
out the sign s^, 2 , we use Eq. d35l l and Table U to show that 
(44) these operators must have four — sT, 2 and two +s\ 2 rotation 



eigenvalues. It immediately follows that 
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Actually, one can recover this result without resorting to nu- 
merics using the following argument. Note first that the quan- 
tum numbers for each single-particle state must be identical 
for the two spin species. Assume that as flux is inserted, no 
single-particle levels cross zero energy, as is typically the case 
in our observations. The quantum numbers for the states be- 
low zero-energy are then conserved under flux insertion. For 
simplicity, let us assume that the half-filled state |0) with no 
added flux carries trivial quantum numbers (which is by no 
means essential). This implies that if for each spin the lower 
zero-mode {i.e., the one pushed downward in energy due to 
finite-size effects) has eigenvalue e za,r ^ 2 under rotation, then 
all other negative-energy states must have eigenvalue e -4 "^/ 2 . 
Denote the upper zero-mode eigenvalue for each spin by 
gVTr/2 One can then easily show that under rotation, the op- 
erators $!| , must have one trivial eigenvalue, one eigenvalue 
e 2i(/3 7r/2 -a 7r/2 ) ) and four eigenvalues e^-/ 2 ""-/ 2 ). The only 

2 > — — 1, which yields 



consistent possibility is for e 4 ^"/ 2 " 



r/2 

Deducing the sign s x is more delicate. To this end we 
consider the composite operation R X T X , which is convenient 
since it does not change the sign of the flux inserted. This 
combination does, however, require a particle-hole transfor- 
mation, so we can not simply read off the eigenvalues of the 
half-filled states from numerics as we did for the rotations. An 
argument similar to the one raised in the previous paragraph 
does nevertheless allow us to make progress. As before, we 
consider a finite-size system where R X T X is a well-defined 
symmetry. A system with periodic boundary conditions along 
the x-direction and hard-wall along the y-direction is partic- 
ularly convenient since one can then insert 2n flux without 
any difficulty. To make the eigenvalues well-defined here, we 
must imagine this flux being inserted slowly so that we can 
monitor the wavefunction continuously during the evolution. 
Assuming no zero-energy level crossings (this has been ver- 
ified in most cases; see below), then there must be at least 
one half-filled state with two zero-modes filled that carries the 
same quantum numbers as the original half-filled ground state 
before the flux insertion. In particular, both states must be 
spin singlets. Now, using Eqs. ( f33l > and d34b along with Table 
U one can readily show that the spin singlet operators $ \ 2 3 
all have eigenvalue s x under R X T X . So we conclude that 



1. 



(46) 



Although we have now fully determined the transformation 
properties of the flux insertion operators $J , it will be use- 
ful to specialize to the 7r-flux state and deduce the sign sq 
that appears under the particle-hole transformation Cq- For 
this purpose we consider the combination T x Cq, which is a 
simple translation whose eigenvalues are easy to determine 
numerically. As above, we consider an L x by L y system with 
periodic boundary conditions along the ^-direction and hard- 
wall along the (/-direction, and choose the Landau gauge for 



all flux configurations so that T x Cq can be realized without 
a gauge transformation. Flux insertions are placed uniformly 
over several consecutive rows midway between the hard walls. 
We restrict ourselves to the case where L x /2 is odd, since the 
"zero-modes" that appear in the presence of 2tt flux can be 
unambiguously identified for such systems. As in our anal- 
ysis of rotations, we examine the six ways of filling the two 
zero-modes, and find numerically that there are four — 1 and 
two +1 eigenvalues under T x Cq for the operators $J . Using 
Eqs. (l33l and (l37b and Table U one can also deduce from our 
earlier results that these operators must have four sq and two 
— sg eigenvalues under T x Cq, implying that 



so 



1. 



(47) 



Note that we have confirmed here that typically there are 
indeed no zero-energy level crossings during flux insertion. 
Moreover, the sign sg can be recovered without numerics us- 
ing the same logic as we outlined for rotations, though we will 
not repeat the argument here. 

The transformation properties for the flux-insertion opera- 
tors $J under all symmetries are summarized in Table [Til 



C. Definition of monopole operators 

We will now define the monopole operators as follows, 

(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
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= 4+-^ 
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rt 


= *L + *6 



Ml = $L+*4,- 



We have organized these "ladder" operators such that the 
monopoles add the same quantum numbers when acting on 
ground states within the q = 0, ±1 monopole charge sectors. 
For instance, A/| adds S z = 1 by filling two spin-up zero- 
modes when acting on |0) and by annihilating two spin-down 
zero-modes when acting on D^ _|— ). Furthermore, these op- 
erators have been defined so that they transform into one an- 
other under the emergent SU(4) symmetry enjoyed by the crit- 
ical theory 20 , implying that all six have the same scaling di- 
mension. Thus the various competing orders captured by the 
monopoles are unified, just as is the case for those encoded in 
the spinon bilinears whose correlations are enhanced by gauge 
fluctuations 20 . Again, this constitutes a highly nontrivial, and 
in principle verifiable, experimental prediction which we will 
elucidate further below. 

Before exploring the competing orders, we note that there 
is another important set of related operators that one should 
consider, which are the following composites involving the 
monopole charge operator Q, 



M] = {M},Q}. 



(54) 



TABLE II: Transformation properties of the flux-insertion operators 
$J . The gauge transformation Co applies only in the 7r-fiux state. 
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[Such operators effectively send $., 



-$,•_ inEqs. gg) 



through (I531 I.1 Our analysis thus far does not enable us to 
distinguish which of these two sets of operators dominates at 
the staggered-flux fixed point. The following argument, how- 
ever, suggests that both sets have the same scaling dimension. 
Consider the current J p = -kz€ iivp F V p, where F vp is the field- 
strength tensor. The monopole charge operator is given by an 
integral over J°: 



Q 



dxdyJ , 



(55) 



which clearly yields an integer q if there is 2irq flux present. 
To all orders in 1/N, J p scales like an inverse length 
squared, 18 implying that Q has zero scaling dimension. Typ- 
ically knowing the scaling dimension of two operators is not 
sufficient to determine the scaling dimension of the compos- 
ite. However, since Q is not a local operator, but rather an 
integral of a charge density, the scaling dimensions for the 
composites MjQ are additive. Thus the scaling dimensions 
for Mj and M.j should be equal. 



IV. COMPETING ORDERS ENCODED IN MONOPOLES 

Now that we have all transformation properties for the flux- 
insertion operators <f>J , we can finally deduce the quantum 
numbers of the six monopole operators defined in Eqs. (l48b 
through (l53b and explore the competing orders encoded in this 
sector of the theory. To this end, we will examine in detail 
the quantum numbers carried by the 12 Hermitian operators 
AfT+Mj wdi(Mj— Mj). These are summarized in TableHTTl 
which is the main result of this paper. (The quantum numbers 
carried by the Hermitian operators constructed from A4j can 
be trivially obtained from these, and we will only comment on 
such operators briefly at the end.) 

In contrast to the monopole scaling dimensions, the am- 
plitudes for their correlations are non-universal and will only 
be related where required by symmetry. We can gain some 
intuition for which operators have the dominant amplitudes, 
at least for weak staggering of the flux, by examining their 
quantum numbers under the particle-hole gauge transforma- 
tion Cq in the 7r-flux ansatz. Those which are even under 
this operation will survive projection into the physical Hilbert 
space, and are thus expected to have the largest amplitudes in 
the staggered-flux case as well. Those which are odd vanish 



upon projection and should have suppressed amplitudes. In 
passing we note that a similar analysis may provide useful, 
though non-universal, information for the flux-conserving op- 
erators as well. The first six Hermitian monopole operators 
listed in Table [Til] are expected to have dominant amplitudes 
by the above logic, while the latter six should be suppressed. 
We proceed now to discuss the results, comparing with previ- 
ous results for the well-studied monopole-free sector 20 where 
appropriate. 

The first operator in Table Hill interestingly, is a singlet that 
carries no nontrivial quantum numbers, and thus constitutes an 
allowed perturbation to the Hamiltonian; we discuss possible 
implications of this in the next section. Note that there is no 
symmetry-equivalent operator in the set of fermionic spinon 
bilinears, all of which carry nontrivial quantum numbers 20 . 
As an aside we comment that naively it may appear, given 
our quantum-number-conservation argument employed ear- 
lier, that having one singlet monopole operator carrying no 
quantum numbers is generic. We stress that this is not the 
case. We applied this argument in different geometries, which 
were designed so that the symmetry under consideration was 
realized in a particularly simple way. Within each geome- 
try, there must be one singlet flux insertion which transforms 
trivially as claimed. But there are three such singlet opera- 
tors, so the same one need not transform trivially in all cases. 
Indeed, similar arguments applied to monopoles on the trian- 
gular lattice yield no such operators carrying trivial quantum 
numbers.— 

Remarkably, the next five operators encode perhaps the 
most natural phases for the square-lattice antiferromagnet — 
valence bond solid (VBS) and Neel orders. We find it quite 
encouraging that these appear as the dominant nearby orders 
in our analysis. Both VBS and Neel fluctuations are also cap- 
tured by enhanced fermion bilinears, which are labeled N c ' 
and N \, respectively, in Ref. 2fJ. It is intriguing to note that 
a recent study that neglected monopoles but took into account 
short-range fermion interactions found that the staggered-flux 
spin liquid may be unstable towards an SO(5)-symmetric fixed 
point, at which Neel and VBS correlations were unified.— In 
light of our results, it would be interesting to revisit that work 
with the inclusion of monopoles, which for the physical value 
N = 4 may also play an important role. 

The remaining six operators in the table are expected to 
have suppressed amplitudes compared to the operators dis- 
cussed above. The first of these transforms microscopically 
like 

M x + h.c. ~ {-iy*+rv [ Sq . (s 6 x S c ) - S b • (S c x S«j) 
+ S c • (S„ x S„) - S d ■ (S„ x S b )], (56) 

where S a = S r _y, Sb = S r+ x, S c = S r+ y, and Sd = 
S r _x. This operator represents the (ir, ir) component of the 
scalar spin chirality. Apart from the finite momentum carried, 
Mi + h.c. carries the same quantum numbers as the enhanced 
fermion bilinear denoted M in Ref. 20 that when added to 
the Hamiltonian drives the system into the Kalmeyer-Laughlin 
spin liquid 31,32 which breaks time-reversal and reflection sym- 
metry. 



TABLE III: Quantum numbers carried by Hermitian monopole operators constructed from M denned in Eqs. J48t through J53l l. In columns 
3 through 5, we provide the eigenvalue if the operator is diagonal under the corresponding symmetry; otherwise the operator into which it 
transforms is given. While all 12 operators have the same scaling dimension, the first six are expected to have the dominant amplitudes based 
on proximity to the 7r-flux state. 
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-1 


-1 


1 


Triplet 


uniform vector spin chirality 



The next two singlet operators in the table transform like 
the following microscopic spin operators, 

M 2 + h.c. ~ (-l) r *[Si ■ (S 2 x S 3 ) - S 2 • (S 3 x S 4 ) 

+ S 3 • (S 4 x Si) - S 4 • (Si x S 2 )] (57) 

Ah + h.c. ~ -(-l) r »[Sr(S 2 xS 3 )-S 2 .(S 3 xS 4 ) 

+ S 3 • (S 4 x Si) - S 4 • (Si xS 2 )], (58) 

where we have used abbreviated notation with Si = S r , 
S 2 = S r+ x, S 3 = S r +x+^, and S4 = S r+ y. These monopole 
operators are closely related to an enhanced fermion bilinear, 
dubbed iV<2, in Ref.|2Cj that transforms like 

Nl ~ Si • (S 2 x S 3 ) - S 2 • (S 3 x S 4 ) 

+ S 3 • (S 4 x SO - S 4 • (S x x S 2 ). (59) 

Furthermore, Ref. |2fJ observed that Njj also possesses the 
same symmetry as the (it, it) component of the skyrmion den- 
sity ps, 



Ps = -7-n- (d x ii x d y n), 

47T 



(60) 



where n is a unit vector encoding slow variations in the Neel 
order parameter. Consequently, A/2,3 + h.c. are symmetry 
equivalent to the (0, ii) and (tt, 0) components of the skyrmion 
density. 

Finally, the last three triplets in the table transform like 
components of the spin operator 



Si x S 3 - S 2 x S 4 



(61) 



Thus, these operators represent the uniform part of the vec- 
tor spin chirality. Enhanced fermion bilinears (N^{ in Ref. 
uOl) also represent vector spin chirality fluctuations, though at 
momenta (0, tt) and (tt, 0). 

What about Hermitian operators constructed from the com- 
posites Aij ? Their quantum numbers can be easily deduced 
from those listed in Table [III] by noting that the monopole 



charge is odd under translations, reflection, and Cq (in the 
7r-flux state), but even under rotations and time-reversal. Con- 
sequently, Hermitian Aij operators have relative momentum 
(tt, tt) and opposite parity under reflection compared with the 
corresponding Mj operators. One can repeat the analysis 
given above for the latter, but we choose not to do so here. 



V. DISCUSSION 

In this paper we have attempted to help resolve an outstand- 
ing issue in the study of algebraic spin liquids — namely, the 
quantum numbers carried by monopole operators — by con- 
sidering the well-studied case of the staggered-flux state. Our 
study builds on previous work 22,23 in the slightly different con- 
text of "algebraic vortex liquids", and can be generalized to 
other settings as well. Essentially, our analysis was predi- 
cated on the assumption that the leading monopole quantum 
numbers can be deduced from the symmetry properties of the 
mean-field ground states with and without a flux insertion, 
with no additional Berry phase effects. While we believe this 
is reasonable, and find the end results to be quite natural, such 
issues can be delicate since we are dealing with a gapless state. 
Thus, we encourage further scrutiny of the conclusions reach 
in this paper. Projected wavefunction studies of the type de- 
scribed in Ref. [33j provide one distinct approach which can 
shed further light on the problem and may help to support our 
findings 34 . A more dynamical treatment of monopoles, how- 
ever, may ultimately be required. 

Assuming we have succeeded in finding the quantum num- 
bers of the leading monopole operators, one issue is worth dis- 
cussing further. Specifically, our analysis showed that there 
is one Hermitian monopole operator which carries no quan- 
tum numbers and thus represents a symmetry-allowed per- 
turbation to the Hamiltonian. An important issue is whether 
this perturbation destabilizes the staggered-flux state for the 
physical number of fermion flavors, which is N = 4. The 
single-monopole scaling dimension computed in the large- N 



limit in Ref.Qjis A m « 0.2657V. Extrapolating to N = 4 
yields A m w 1.06, substantially lower than 3, suggesting 
that the symmetry-allowed monopole operator may constitute 
a relevant perturbation. However, caution is warranted here 
(more so than usual in such extrapolations), since the sub- 
leading correction to the scaling dimension is generically an 
TV-independent, possibly 0(1) number. Given the obvious 
importance of this question for the high-T c problem, further 
studies of these scaling dimensions are certainly worthwhile. 
And if the operator turns out to be relevant, what are the prop- 
erties of the phase to which the system eventually flows? An 
interesting possibility is that the system may flow off to a dis- 
tinct spin liquid state, but it is also possible that dangerously 



irrelevant operators lead to broken symmetries. This question 
is left for future work. 
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